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Abstract 

We study particle decay as the origin of dark radiation. After elaborating general 
properties and useful parametrisations we provide model-independent and easy-to- 
use constraints from nucleosynthesis, the cosmic microwave background and struc- 
ture formation. Bounds on branching ratios and mass hierarchies depend in a unique 
way on the time of decay. We demonstrate their power to exclude well-motivated 
scenarios taking the example of the lightest ordinary sparticle decaying into the 
gravitino. We point out signatures and opportunities in cosmological observations 
and structure formation. For example, if there are two dark decay modes, dark ra- 
diation and the observed dark matter with adjustable free-streaming can originate 
from the same decay, solving small-scale problems of structure formation. Hot dark 
matter mimicking a neutrino mass scale as deduced from cosmological observations 
can arise and possibly be distinguished after a discovery. Our results can be used 
as a guideline for model building. 



1 Introduction 



New cosmological probes measure the amount of radiation in the Universe at different 
epochs of its evolution with a crucial increase in precision. One strength of the standard 
cosmological model, which is based upon general relativity and the Standard Model 
(SM) of particle physics amended by "invisible" components known as dark matter and 
dark energy, is the precise prediction of the amount of radiation. Within the first 20 
minutes light nuclei like helium were formed during the process of big bang nucleosyn- 
thesis (BBN) as observed today. At such early times the Universe was dominated by 
radiation. Since nucleosynthesis depends on the expansion rate, BBN is sensitive to 
the amount of radiation. There is still a controversy between different groups, some 
favouring the prediction and others an increased amount [1, 2, 3]. Since the main uncer- 
tainty stems from the determination of the relic helium abundance from observations, 
we can expect improvements in the foreseeable future. Observations of the cosmic mi- 
crowave background (CMB) constrain the amount of radiation in the Universe in an 
epoch lasting from some thousand years till photons decouple 10 5 years later. Since 
the first determination of the radiation content of the Universe from the CMB roughly 
ten years ago [4] and also in current measurements by the South Pole Telescope [5] and 
the Atacama Cosmology Telescope [6], mean values are larger than the prediction. The 
observed suppression of the CMB power spectrum at larger multipoles would be due to 
increased Silk damping [7]. Non-Gaussianities could provide further insights [8]. Ad- 
ditional radiation also eliminates tension between cosmological data and measurements 
of today's expansion rate [9, 10]. Increased mean values are found as well in extended 
analyses including additional cosmological data [11, 12, 13, 14, 15, 16, 17, 18], although 
the results for the statistical significance of this increase vary. More importantly, due to 
the increase in precision [19, 20] the Planck satellite, which finished data taking already, 
could turn these hints into a 3-er to 5-cr discovery, if current mean values are accurate. 
Our understanding of the third component in the Universe -besides matter and vacuum 
energy- would be proven incomplete, too. Since such additional, "invisible" radiation 
cannot arise from the SM and its nature is unknown, it has been dubbed dark radiation. 

In this work we study particle decay as the origin of dark radiation. Additional 
radiation has been studied mainly in connection with the possible existence of additional 
neutrino species beyond the known three [21, 22, 23, 24] or other relativistic species [25, 
26, 27]. The amount of additional radiation is then generically discrete. The species 
are relativistic during BBN and may or may not still be relativistic around photon 
decoupling. Hence, they always lead to an increase in radiation at BBN and generically 
the same increase during CMB times. This is an appealing prediction, in particular, 
because past observations have been too imprecise to find a difference between BBN 
and CMB determinations. In contrast, dark radiation from particle decay can originate 
before [28, 29, 30, 31, 32, 33], during or after BBN [34, 35, 36, 37, 38, 39, 40, 41] and can 
become non-relativistic before, during or after photon decoupling. Any observed value 
can be explained. These are qualitative differences to the case of a relativistic species 
that become distinguishable given the new observational precision. Even though current 
observations do not allow for any conclusion [7, 42], the most interesting finding would 



2 



be an increase in radiation after BBN. This would exclude additional relativistic species 
leaving late particle decay as the most attractive possibility. 1 Interestingly, the time of 
decay can be probed in the inflationary gravitational wave background [44]. 

We exploit the fact that the energy density of the decaying particle is fixed by the 
observed amount of dark radiation to determine model-independent upper bounds on 
several branching ratios of the decaying particle from BBN, spectral distortions in the 
CMB and the ionisation history of the Universe. Thereby, we point out an opportunity to 
solve the cosmological lithium problems [45] and the discovery potential of a future CMB 
polarimeter for the considered decay. A decay before BBN could mimic a cosmology with 
additional relativistic species. More importantly, there is a plethora of new cosmologies. 
We elaborate constraints and opportunities from heavier decay products in structure 
formation. They may just form dark radiation, but they do not need to act as radiation at 
all. If they form the observed dark matter, two of three dark components would originate 
from the same decay. If they are not cold, their free-streaming might resolve the missing 
satellites problem [46, 47]. While lighter decay products act as dark radiation, heavier 
ones might mimic the neutrino mass scale as deduced from cosmological observations. 
Cluster abundances seem to favour additional radiation together with a finite neutrino 
mass scale [48]. 

In the next section we study the simplest case allowing for exactly one dark decay 
mode. We will use our findings in Sec. 4 to explore which opportunities open up in 
structure formation, if there are two dark decay modes. Sec. 3 is devoted to general 
constraints and opportunities from BBN and CMB. We summarise and conclude in 
Sec. 5. In the appendix we provide an analytic treatment of the exponential decay law 
in an expanding universe. 

2 One dark decay mode 

In this section we study the origin of dark radiation from (one) two-body decay of a 
non-relativistic particle. Such a decay is drawn in the upper right corner of Fig. 1, 
where we indicate some nomenclature. The decaying particle (mother) decays into non- 
identical particles, where one is heavier (heavier daughter) and the other one lighter 
(lighter daughter). In Fig. 1 we illustrate the behaviour of certain energy densities 
in a cosmology with dark radiation from particle decay. First of all, the plethora of 
possible cosmologies cannot be shown in only one figure. So we illustrate a few typical 
and interesting scenarios. At some early time corresponding to a small scale factor the 
Universe is dominated by radiation and the mother is relativistic, so that their energy 
densities scale equally. At some time the mother becomes non-relativistic and from then 
on its energy density p tx a -3 grows relative to the radiation energy density p ia _d oc a~ 4 
as the Universe expands, where a denotes the growing scale factor. Thus stable or very 
long-lived matter generically comes to dominate the Universe. Actually, this fact gives 
rise to various cosmological problems with the well-known gravitino problem as prime 
example. These problems may turn out as fortunes, if they give rise to the desired 

One alternative is the late annihilation of additional species [43] . 
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Figure 1: Behaviour of comoving energy densities pa? in an expanding universe with 
dark radiation from particle decay. The fulhlogarithmic figure is illustrative and not 
exact. Upper right corner: Nomenclature for the considered two-body decay. 



dark radiation as in [36]. Since only the relative behaviour matters, comoving energy 
densities pa 3 are drawn. When the mother decays, its energy density is converted into 
the energy density of its relativistic daughters. The radiation content of the Universe is 
increased. Of course, there are various possible production mechanisms for the mother 
in the early universe. It is crucial only that its energy density is within a certain range 
at its decay, see below. The energy density of the daughters scales as radiation till they 
become non-relativistic. They may still be relativistic today or, especially, the heavier 
one may have become non-relativistic earlier and thus could possibly form the observed 
dark matter or some hot dark matter component. For comparison, the energy densities 
of radiation and matter in standard cosmology are shown as grey dashed curves. 

If the decaying particle (mother) decays with some effective strength y into its decay 
products (daughters), an effective decay width might be given as T ~ y 2 m/(167r), where 
m denotes the mass of the mother. For the following it is assumed that the branching 
ratio of this decay is close to one. Actually, we find in Sec. 3 that in all cases we are aware 
of the branching ratio into dark components is constrained to be very close to one at 
times later than £bbn ~ 0.1 s. Such branching ratios are common in dark matter models, 
because usually some symmetry is invoked to stabilise the dark matter candidate. In 
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the on-shell tree approximation some heavier particle sharing the symmetry then has to 
decay into the dark matter candidate. Such branching ratios may also -or in addition- 
be enforced by the mass spectrum allowing for only one decay channel kinematically. 

A prime example for a dark matter stabilising symmetry is -R-parity in supersym- 
metric models, which also naturally comprise extremely long-lived particles, if combined 
with gravity. Think about the gravitino decaying into axino and axion with effective 
y 2 ~ "i|/ 2 / where m 3 / 2 denotes the gravitino mass and M p \ the reduced Planck 
mass. This decay naturally leads to the emergence of dark radiation way after BBN but 
before photon decoupling [36]. The mother (= gravitino) would decay into a fermionic 
axino and an axion scalar. Another example is the decay of the lightest ordinary su- 
persymmetric particle (LOSP) into its superpartner and the gravitino with effective 
y 2 ~ m fosp/(3-^pi m 3/2)' wnere mi osp denotes the LOSP mass. The scenario reminds 
of the sWIMP mechanism, where these kind of decays were considered to produce the 
observed dark matter. In the case of a neutralino LOSP the mother were a fermion de- 
caying into the fermionic gravitino and a gauge boson. In the case of a sneutrino LOSP 
the mother were a scalar and both daughters (= neutrino and gravitino) fermions. We 
mention a third example. In higher dimensional theories the superpartner of a modulus 
field, i.e., a modulino might decay with effective y 2 ~ ^ m ^/ where denotes the 
modulino mass and A an effective coupling. The modulino might decay into a sneutrino- 
neutrino or axino-axion pair and so on. Note that there are various combinations of 
spins. 

2.1 Basics 

We introduce useful parameters, determine general properties and derive basic equations. 

Kinematics In the rest frame of the decaying particle the decay products of a two- 
body decay have in full generality momenta with opposite direction and same absolute 
value. It is 

\pi\ = \P2\ = y~ (( m2 ~~ ( mi + m 2) 2 ){m 2 - (mi - m 2 ) 2 )) 2 , (1) 

if m denotes the rest mass of the decaying particle and subscripts 1,2 label the two decay 
products. We choose subscripts such that mi < mi. We find it useful to define 

s ^rn-m 2= rn L _ 1>Q 
ni2 rri2 

as measure of the mass hierarchy between mother and the heavier daughter; or their mass 
degeneracy for S < 1. Negative 5 does not exist, because the decay were kinematically 
forbidden in that case. The mass of the heavier daughter can be written as 

mi = (S + iy 1 m. (3) 

If the decay shall produce dark radiation, the energy of the lighter decay product, 
E\ = y/\p{ \ 2 + m 2 , must be dominated by its kinetic energy. For mi <C \p±\ its en- 
ergy E\ ~ \pt\- In the limit mi <C mi, which is equal to mi/m <C (5 + 1) — 1 , the general 
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momentum (1) simplifies as 



m(<5 + l) 2 -l m (. 1 \ ... 

m \T m2 W = 2 (5 + 1)2 = y ^ - (^2 J • W 

For the case mi = mi see Sec. 4. 

On T" r An initial particle momentum p- in i from a decay at temperature decreases 
due to the expansion of the Universe. Here and in the following, temperatures T refer 
to the corresponding photon temperature T 7 at the considered time, e.g., = T 7 (r) = 
T(r). The momentum at temperature T is 

i 



p(T) = Pini— = Pini— ( -j- , (5) 



where a denotes the scale factor of the Universe and the second equality is due to the 
conservation of comoving entropy. As usual, g* s = g* s {T) is the effective number of 
degrees of freedom in the entropy density of the Universe s = (2tt 2 /45)<7* S T 3 . The 
superscript on g* s indicates here and in the following at which temperature or time 
g* s is evaluated, g* s = g* s (Td)- The same holds for subscripts on a. We define the 
temperature T nr , when a particle becomes non-relativistic, by 

p(T nr )=m, (6) 

where p = \~p*\ denotes the particle's momentum and m its mass. For a particle species 
following the distribution P we consider the mean momentum to determine whether the 
species is relativistic or non-relativistic. 

In our case p m \ is given by (4) and the mass of the heavier daughter by (3). Therefore, 
the condition 

p 2 (T?) = m 2 (7) 

yields 

with g™ = ^(T"). The correction factor fj, = c _1 r[c _1 ] takes into account the expo- 
nential decay law in an expanding universe, a oc t l / c , compared to the sudden decay 
approximation. It is derived in Appendix A. It is p = n(P) — 0.886, if the decay occurs 
during radiation domination and ji ~ 0.902, if the decay occurs during matter domina- 
tion. Throughout this work we will often argue under the assumption of a sudden decay, 
because this simplifies the discussion and reveals key points. One example for this is 
whenever we consider some notion "at decay". We will take into account corrections 
from the exponential decay law in the final equations by correction factors, which rep- 
resent good approximations for times t > 3r or > 4r , cf. Appendix A. Often times of 
interest are indeed much later than the time of decay. 
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Energy densities The "non-dark" radiation energy density of the Universe, i.e., the 
energy density of the thermal bath in the Universe, is given by 

J2 



Prad 



7T~ 

30 



(9) 



where g* = g*(T) denotes the effective number of relativistic degrees of freedom in the 
bath. Bounds on the total radiation energy density p\°^ d exist from processes around 
and after e + e~-annihilation, so for cosmic temperatures around and smaller than the 
e + e~ -annihilation temperature T e + e - ~ m e ~ 0.5 MeV. They are usually given in terms 
of the effective number of neutrino species N e s defined by 



tot 
Prad 



1 + ^Veff 




(10) 



where the radiation energy density is given as a sum of the energy density in photons 
p-y = (vr 2 /15)T 4 , the energy density in SM neutrinos with A^m = 3.046 [49] and T u /T^ = 
(4/11) 1 / 3 and any departure from the standard scenario parametrised as a summand 
in N c g = + AN e g. The small correction in is due to incomplete neutrino 

decoupling at e + e~-annihilation. We denote temperatures before neutrinos become non- 
relativistic at T° r and lower than T e + e - by lj ow . Comparing (9) and (10) we see that 
9*{Tiow) — 3. 384. 2 In this temperature range the energy density in dark radiation p^ T 
can be written as 

p dT = AN cS x p lv = AN cS x I ( ^ J p 7 , (11) 



if pi u denotes the energy density of one SM neutrino species with thermal spectrum. It 
follows that 

Air 



0.1342 x AiV cff 



(12) 



low 



for T° r < T < T e + e -. We see that even for the 5-a exclusion, AA^ ax 



5.265, of the 

combined analysis in [6], there would be less dark than SM radiation. Here and in the 
following a vertical line with subscripts like | low indicates at what time (or temperature) 
the corresponding term is evaluated. Towards higher temperatures the bath energy 
density p ra d oc g*T 4 , while the one in dark radiation scales as pd r oc gfi 3 T 4 , if dark 
radiation is decoupled from the bath. Thus, at the time of the decay producing the dark 
radiation, 



Pdr 



Prad 



dec 



Pdr(^dcc) Prad(^l ow) Pd 
Pdr Plow) Prad(?c 



dec ; 



SM 
Prad 



low 



Q d 
^low 



low 



Pdr 



SM 
Prad 



low 



and inserting g*(7i ow ) and g*s(T\ ow ) 



y*. 



0.5440 



2(1 + A^28/88) ~ 3.938 finally 2 

(<d) 4/3 Pdr 



SM 
Prad 



(13) 



(14) 



low 



' The differences in 

jSM 
Veil • 



and g, s to the often used values in the literature are due to the small correction 
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which is valid at any decay temperature. Due to the different scaling behaviour the dark 
radiation component could have dominated the Universe at decay, but only for decays 
with = g* s 3> and for extreme values of AN e s- At intermediate temperatures T 
we consider the different scaling of dark and SM radiation. It is 

o 4/3 

Pdl (T) = 0.0730AiV eff ^p rad (r) . (15) 
9* 

This is a useful parametrisation of the dark radiation energy density. 

The desired amount of dark radiation determines the energy density of the decaying 
particle p = n m at its decay. In a two-body decay with branching ratio one the number 
densities of the decay products are fixed to be equal to the number density of the decaying 
particle, n = n\ = 112- As radiation the energy of a particle can be approximated by its 
momentum, E ~ p, and the heavier daughter may act or may not act as dark radiation 
at the times of observation. Then the energy density of dark radiation at decay reads 
Pdrldec — 9^ Sn Ei ^ g^npi, where p± is determined from the kinematics (4) and g^ s 
counts the number of dark radiation components at the time of observation. It is gd T = 2 
or 1, if the heavier daughter particle became non-relativistic before the time probed by 
observations. We define a conversion factor / by 

Pdr = f x p (16) 

so that 



gg s (S + l) 2 -l 
2 (6 + 1 

The factor p takes into account the exponential decay law in an expanding Universe, see 
Appendix A. Altogether, we obtain 

0.146 AiV eff (<4) 4 / 3 (5 + V 2 



P 



Prad 



! P,r 



dec Arad 



dec P 9t* 9* [5 + 1? -I' (18) 



We see that the decaying particle could be required to dominate the energy density of the 
Universe for a short time prior its decay to explain extreme values of AN c s. However, 
this is not to be expected and, especially, at late times improbable. The dilution relative 
to p ra d is only oc gl^ 3 , so that a large number g* S> g^ M would be necessary to change 
this. From (18) follows today's energy density of the decaying particle, if it were not 
decayed, in units of today's critical energy density p c as 

™ 2 = %^(|) 8 f, (19) 



where the dilution till today, that would have happened, is considered. In this work we 
will heavily use the fact that the energy density of the decaying particle is fixed by the 
amount of dark radiation and thus by observations independent of an underlying particle 
physics model. The other way around we can single out AN e g in (18) resulting in 

AN cS - 9.13/,— (5 + i)2 ( ^ )1/3 , (20) 
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where we introduced the particle yield Y = n/s. 

If the decay occurs before BBN, r <C £bbn, the decaying particle may dominate the 
Universe prior its decay. In this case AN e g is set by the relative branching into radiation 
formed by SM particles and dark radiation. This relative branching is just found in (15), 
i.e. 0.073 AN cS gl /3 . 



2.2 5-t plane 

The mass hierarchy between decaying particle and decay products is constrained from 
several cosmological considerations. First of all, at no time there is an upper bound on 
8, because the energy density of the decaying particle can, in principle, be adjusted such 
that all its energy transferred to radiation at its decay accounts for the desired increase 
in AiV e g. This might be different in a concrete particle physics model, where 5 and/or 
p turn out to be given. More importantly, there are lower bounds on 5 from cosmology. 
The underlying considerations are qualitatively different, if the decay occurs before or 
after matter-radiation equality at t eq . In any case they depend on the time of decay r, 
while they are independent of the underlying particle physics model. Our assumptions 
on the time and temperature of equality are outlined at the beginning of Sec. 2.3. 

Before i eq there is a lower bound on 5 from the requirement that no daughter particle 
may come to dominate the Universe before t eq . This non-dominance requirement can be 
expressed as 

VL2h 2 = frmax^dm^ 2 (21) 

with 6 max < 1. When the mother decays (not too close before t eq ) also the heavier 
daughter particle must be emitted as radiation. Otherwise, it would dominate already 
shortly after being emitted. Thus its energy density is as large as the energy density of 
the lighter daughter, p2\dec = Pi Idee = (5dr >S ) _1 / 9 dr(2d). After emission its energy density 
scales as radiation oc a -4 till it becomes non-relativistic at . If T 2 nv < T eq 43- > t eq , 
P2 surely never dominates. From becoming non-relativistic on it scales as matter oc a~ 3 . 
Thus 



PdrCTd), obs . ^T* (gl 



n rurv-uy / bs\-l 'l I / 00 n 

n2= — {9 ^ ) ^rUJ ' (22) 

where we anticipated = g® s . Inserting p^ from (15) the inequality (21) becomes an 
upper bound on T^ r 

T 2 nr < 6.938 Tcq AN^b max ^^gt s (23) 

where we can see that for the bound only g^ s = 1 is sensible. Since T^ r <C T e + e -, it is 
justified to set g™ = g®. We find the corresponding cosmic time as 

(^ r ) min = 0.0254 t cq b^ x AN* s (^§) 2 = 4.61 x 10 10 s6 m L AiV e 2 ff Q^) 2 • (24) 

For times later than (iJJOmin the assumption of relativistic emission is no longer neces- 
sarily fulfilled. We note that from (23) is likely later than the time when the first 



9 



5 m in(r) from different requirements 




T[S] 



Fi gure 2: <5 m i n -T-plane exploiting (25). Values above the corresponding line are consid- 
ered to be allowed. The thick solid curve corresponds to the hot dark matter constraint, 
cf. Sec. 2.3.2, with AN e g = 1. Thin dashed curves below (AiV c fj = 0.52) and above 
(AN c ff = 5.265) show the dependence of this bound on the produced amount of dark ra- 
diation. The thin solid curve corresponds to the non-domination constraint. The dotted 
curve (at the lower edge of the thick solid one) gives 5 such that the heavier daughter 
becomes non-relativistic at photon decoupling i^ pl and the dash-dotted one such that 
this happens today. These three curves are for AN c q = 1. At the upper edge of the 
thick solid curve mean values of [48] originate from the decay for massless neutrinos. 
Various important and suggestive times are highlighted by vertical dashed lines: onset 
tbbn and end of BBN i^bn' the earliest possible time for the heavier daughter to become 
non-relativistic (i^Omin from (24), re-entry of first observable modes in the CMB t™£ 
and matter-radiation equality i eq . For decays during BBN the increase in N e s deter- 
mined from BBN is smaller than the corresponding increase measured in the CMB. The 
relative difference depends on the time of decay as quantified in [38]. In all figures we 
take into account the evolution of g* and g* s . That the curves are smooth around t e + e - 
shows that dependencies have cancelled. Within the horizontal dashed lines the mass 
hierarchy or degeneracy is within an order of magnitude. 
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observable modes of the CMB enter the horizon, which sets t™£ ~ 5.2 x 10 10 s [35]. 
Indeed, this is for sure taking into account constraints from structure formation requir- 
ing 6 max < 1, cf. Sec. 2.3. The heavier daughter is restricted to become non-relativistic 
during CMB times or later, which might leave observable consequences due to the cor- 
responding change in the expansion rate. 

Inserting from (8) in (23) we obtain an implicit lower bound on 5 



i 



If 6 3> 1, the lhs of (25) reduces as 

1 (<5 + l) 2 - 1 

lim - ' = 1. (26) 



<5->oo 5 5+1 



Then (25) becomes practically 



. 0.2883 /t cq V A A7 /0.1246\ f g°Y f gtY 



where we used the time-temperature relation in a radiation dominated universe to replace 
T d /T eq = (<7?7<7* d ) 1/4 (WT) 1/2 ^d g c * = fi. 

After i e q a relativistically emitted, non-dominating particle becomes even more 
subdominant as the Universe expands. However, there is a lower bound on 5, if we 
require some significant AiV c fj > 0, because the maximally allowed energy density for 
the heavier daughter particle is the dark matter energy density as (21) with the crucial 
difference that it is emitted non-relativistically, here. When b ~ 1, dark matter would 
decay and convert a small amount of its energy into radiation. This is qualitatively 
different for r < (i-fOmin- In the intermediate regime, (ifOmin < t i$ t eq , we expect 
observable consequences in the CMB. 

The energy density of a non-relativistic, p2 < m2, species can be approximated as 
p2 = n2-&2 — Again exploiting n = n<i = m and replacing mi by (3) we find 

nro ^ = iy i f -i pdr{Td) , (28 ) 

+ 1 

where we used that p = nm and the definition of the conversion factor (16). Taking 
into account the red-shift till today the energy density of the heavier daughter in units 
of today's critical energy density is given by 



From the requirement (21) we obtain the very same implicit lower bound (25) on 5 for 
decays after (iJfOmin, where only g* s /g® s has to be replaced by one. If S — > 0, the lhs 
of (25) reduces as 

, 1 ((5 + l) 2 -l 

lim — 1 -. '- = 1. 30 

s^o 26 5 + 1 J 
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Then (25) becomes practically 
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. 0.2883 ft cq \ s A Ar /0.1246\ 

where we used the time-temperature relation in a matter dominated universe with con- 
stant g* to replace Td/T eq = (i eq /r)3. Whenever necessary we assume a sudden transi- 
tion from radiation to matter domination in analytic calculations. 

The lower bound on 5 from (25) with AN c g = 1 is depicted in Fig. 2 as thin solid 
curve. The analytic approximations (27) and (31) are overplotted as very thin grey curve 
with a jump and its largest deviation at t eq and a tiny underestimation at very late times. 
A decay before ibbn increases AN e g before BBN. We found that the decay products may 
not become non-relativistic before a time (if r ) m in > T ™h- For decays during BBN 
the increase in iV e g determined from BBN is smaller than the increase measured in the 
CMB depending on the time of decay [38] . Decays after i^n add radiation during CMB 
times. It is usually assumed that the decay products are still relativistic today or at 
least till photon decoupling at t^ pl , while this must not be the case, especially not for 
the heavier daughter. The figure ranges beyond t^ pl to times as late as 10 17 s. This is 
for completeness. We point out that for such late decays the meaning of the curve is 
maybe far from clear in this and other following graphs. This is due to the currently 
limited understanding. Nevertheless, first pieces of information are provided that might 
motivate focused studies of such cosmologies. 

Finally, let us give the following comment. We could have started from the general 
requirement (21) for any time and would have obtained the same results. We would like 
repeat that the cosmology is very different for a decay close to the boundary occurring 
sufficiently before or after (tifWn- After (iifOmin the energy density of the decaying 
particle is of the order of the energy density of the heavier daughter. In some sense 
dark matter decays into today's dark matter particle and dark radiation. Before (t5 r )min 
the energy density of the decaying particle has to be much larger than the dark matter 
energy density to produce significant dark radiation. One could expect a lower bound on 
5 from requiring some significant AiV e g to arise also before (ijfOmin- However, this is not 
the case. The energy density of the decaying particle could in principle be arbitrarily 
large, such that in its decay a tiny amount of its energy converted into dark radiation 
suffices. This is because dark radiation does not thermalise and thus does not dilute pre- 
existing abundances. The additional entropy is negligible as the momenta of the decay 
products are distributed on a small shell in phase space, cf. Appendix A. An arbitrarily 
small 5 seems allowed. Indeed, the non-dominance requirement of (21) were violated. 
For 5s fulfilling (25) the upper bound on the energy density of the decaying particle is 
set by the possible overproduction of dark radiation, so a too large increase in AN c q 
inconsistent with observations. In the scenario under consideration this upper bound 
were set via (18). The lower bound on 5 requiring some minimal AN e g is obtained by 
inserting this value into (25). For times around t eq the situation is more involved than 
outlined, here, because quite a large fraction of the dominating matter component would 
have to decay into radiation. However, we expect such a situation to leave pronounced 
signatures in the CMB. 
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2.3 Structure formation 



The transition from radiation to matter domination is one of the most important events 
in structure formation. The red-shift of matter-radiation equality z cq is related to the 
dark radiation and matter density fi m by (cp. (53) of [50]) 



Naively, one could conclude that AiV e g > thus implies smaller z eq , i.e. a later tran- 
sition. However, this is not observed. Primordial abundances of light nuclei are simply 
affected by the expansion rate during BBN, which is given by the amount of radiation 
and thus N e g. What can be seen in the CMB power spectrum is a suppression at higher 
multipoles (corresponding to smaller scales) compared to the expectation from the stan- 
dard scenario. Constraints on iV e g depend on the cosmological model, which can be 
extended in other ways to suppress small-scale power. For an increase in N e g the sup- 
pression has been identified to be due to increased Silk damping caused by the increased 
expansion rate [7]. Indeed, (32) shows the degeneracy of N e g and Q m h 2 , in particular, 
for WMAP using the first and third acoustic peak to determine z eq . This degeneracy is 
broken by including measurements on smaller scales. The combined data then allows to 
measure iV e ff in addition to z eq , see references in the introduction. The studies do not 
find that z eq varies with AiV e ff. Even though they also do not exclude a smaller z eq , we 
adopt the PDG mean value z eq = 3200 ± 130 [51], which is consistent with the other 
studies. The temperature at equality is thus T eq = (1 + -z eq )To ~ 0.752 eV. We fix the 
time of matter-radiation equality to t eq ~ 1.81 x 10 12 s from the relation between time 
and red-shift in a universe filled with radiation and matter using PDG mean values for 
the time and red-shift of decoupling [51]. 

For a fixed z eq the matter and radiation density are no longer independent. In Fig. 1 
this might be seen comparing the energy density of the Universe including dark radiation 
(thick, solid) with standard cosmology (grey, dashed). The baryon density is regarded 
as robustly measured, Q^h 2 = 0.0226(6) with la uncertainty in the last digit [51]. 
Therefore, the uncertainty in AA^g- turns into an uncertainty in the dark matter density 
as f2 m = fib + ^dm- Following (32) dark matter and dark radiation energy densities are 
linked as 



Omitting this dependence in, for example, Sec. 2.2 can lead to deviations up to few ten 
percent. Thus it is important to take this dependence into account. This is especially 
true when considering more complete cosmologies including, for example, an origin of 
dark matter as in Sec. 4. Particle physics parameters can be significantly affected. An 
example might be the determination of the relic density of the lightest supersymmetric 
particle in supersymmetric models. 

Increasing the dark matter energy density, while keeping the baryon density constant, 
decreases the baryonic matter fraction Qb/^m- This decreases the pressure support on 
matter prior to photon decoupling and, therefore, boosts the growth of structures be- 
low the sound horizon at photon decoupling of about 150 Mpc [52]. A finite neutrino 




(32) 



Q dm h 2 = 0.1067 + 0.0179AA cff . 



(33) 
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mass scale has the opposite effect, because neutrinos then count to Q m , while they do 
not take part in structure formation below their free-streaming scale. Consequently, a 
best-fit might be found having both, additional radiation and a -small, "compensating" - 
dark matter fraction with very large free-streaming scale, see [53, 54] considering sterile 
neutrinos. In this line of thought it is very interesting that AN c g > together with a 
finite neutrino mass scale seems to be favoured by measurements of the abundance of 
galaxy clusters corresponding to ~ 10 Mpc scales [48]. We note that all these consider- 
ations assume a decay before the affected epoch, so for example before i eq to increase 
the amount of radiation at t e q and so on. It is neither clear how an increase in N e s after 
i^ pl might be observed nor what is the effect of the corresponding amount of matter 
being transformed into radiation. Also in this sense some figures range beyond tj P for 
completeness only. 

In the following we shall point out the constraints and opportunities arising from 
the heavier daughter in structure formation. It may constitute just half of the dark 
radiation, but it does not need to act as radiation at all. The opposite extreme case 
would be that the heavier daughter forms the observed dark matter. In between these 
two cases its free-streaming could mimic the finite neutrino mass scale as deduced from 
cosmological observations. 

2.3.1 Minimal free-streaming scale 



(^2)min - ^(^min) 




10 9 10" 10 13 10 15 10' 

r[s] 



Figure 3: Minimal free-streaming scale of the heavier daughter (Ajjmin depending on 
the time of decay. Times are highlighted as in Fig. 2. Horizontal dashed lines indicate 
the galaxy cluster scale A gc ~ 10 Mpc and 1 Mpc as suggestive free-streaming scale of 
warm dark matter, respectively. 
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As collisionless particle the heavier daughter can stream out of overdense regions 
and into underdense regions, in the process smoothing out inhomogeneities. In order to 
take this effect into account properly, one should integrate Boltzmann equations. In the 
following we shall estimate the scale of collisionless damping analytically. The decisive 
quantity is the free-streaming scale 

Af = T V -^dt (34) 



of the heavier daughter, where v(a) denotes scale factor dependent velocities. The situ- 
ation is qualitatively different whether the decay occurs before or after (iJfOmin- 

Before (t^Omin both decay products are emitted relativistically. If a particle be- 
comes non-relativistic after i eq , its free-streaming scale is much larger than the size of a 
protogalaxy. It is said to act as hot dark matter, see Sec. 2.3.2. We will be interested 
in the case used in Sec. 2.2 to determine <5 m j n where the heavier daughter becomes non- 
relativistic before t eq . We approximate the velocity of a particle emitted with relativistic 
momentum by 

.(a)J 1 (35) 
\^ ,ifa>a nr . 

After the emission with v ~ 1 the velocity decreases with the expansion of the Universe 
and becomes smaller than one, when the particle becomes non-relativistic at ttnr* A 
finite, non-relativistic velocity of the decaying particle is safely negligible. Exploiting 
the additivity of integration on intervals we thus find 

A fe « /"*" htt + f " °^dt + r °^dt . (36) 

The first integral corresponds to relativistic free-streaming, the second one to non- 
relativistic free-streaming before t eq and the third to non-relativistic free-streaming dur- 
ing the matter dominated era. We can safely neglect the current vacuum dominated 
phase. 

In order to perform the integration analytically we take g* and g* s to be constant 
for times as late as t nr . As the variation of gli 3 and g\^ with time is weak, we treat 
them like constant factors in the first integral, too. We checked numerically that this 
comes with negligible errors. With the additional approximation of a sudden transition 
between radiation and matter domination at t e q the integrations become straightforward 
exploiting the common a(t) relations. We obtain 

with g^l = <?* s (tnr) and g™ = g*(t nT ). The first term corresponds to relativistic free- 
streaming, the first summand in the brackets of the second term to non-relativistic free- 
streaming during the matter dominated era and the logarithm to non-relativistic free- 
streaming before i eq . Accepting an error of percent we have taken (i e qAo) 1//3 — >• in the 
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last integral. In (37) we can see that for the often discussed case r <C t nr <C t cq relativistic 
free-streaming is negligible and non-relativistic free-streaming before t eq dominates over 
free-streaming after t eq . In this case Af s is roughly set by i nr . In order to investigate 
our case we set t nr = i?> r = ^C^f) m (37). The scale factors at corresponding times 
are given by a d = a (T /T d )(g° s / g* s ) 1/3 and a nr = ao(T /T^)(g° s /g^ s ) 1 / 3 , respectively. 
Since is given by (8) we can express the free-streaming scale of the heavier daughter 
as function of 5 and the time of decay r. We find 



A*(<S,t) 
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(38) 



This analytic approximation applies to any case with r < t nr < i cq <C to- Of course, as 
r — > i nr there is no relativistic free-streaming. 

The minimal free-streaming scale of the heavier daughter (Af^min is indeed inde- 
pendent of the time of decay. It is instead given by the earliest possible time for the 
heavier daughter to become non-relativistic (23), which corresponds to the lowest possi- 
ble initial velocity. It is achieved for minimal 5 given in (25), because for the minimally 
required mass hierarchy the heavier daughter is emitted with its minimally allowed initial 
momentum. We find 

aAk.D « 9 Mp< AA (11 ( 

21 \ 



(A2 )miu i 



+ In 



48 AN. 



o h 2 \ 

1 1 'dm' 1 \ 



eff 



0.1246 J 



(39) 



Consequently it depends on the amount of dark radiation only, except for the relativistic 
free-streaming scale in the first line. It is independent of and g* s at decay. Taking 
AiV c ff = 1 and the limit r <C t nv we obtain (A| s ) m i n ~ 89 Mpc as shown in Fig. 3, which 
should be compared to the galaxy cluster scale A gc ~ 10 Mpc. It is actually roughly a 
factor of nine larger. Thus the heavier daughter from a decay before (t^Ornm is in any 
case way too warm to form the observed dark matter, see below. 

After t e q to act as matter the heavier daughter should be emitted non-relativistically. 
We approximate its velocity by 

Pini Qd 
m,2 a ' 

where p m \ is given by (4). With a oc 

£2/3 

in a matter dominated universe the scale 
factor at decay ad = (r/io) 2 / 3 and the integration in (34) becomes straightforward. The 
free-streaming scale of the heavier daughter for r > t eq reads 



^2(0) 



(40) 
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(41) 
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where, now, 5s are smaller than one in any case. Inserting the minimal 5 from (25) we 
obtain the minimal free-streaming scale of the heavier daughter (A2 S ) m j n for r > t eq . It 
is depicted in Fig. 3. For decays occurring in the small window after (t^Omin 

but before 

t eq the curve in Fig. 3 represents a lower bound, because we did not consider that the 
Universe is radiation dominated in this period. However, around this time, for example, 
the assumption of sudden transitions leads to quite large errors in the estimate. We find 
that a correction were small and thus not important for our purpose. 

Lyman-a forest data constrains the scale of collisionless damping. Constraints range 
between Af s < 0.5 Mpc [55, 56], tighter ones [57] and those more relaxed due to the 
rejection of less reliable data [58, 57, 59]. We assume a constraint on the free-streaming 
scale Af s < 1 Mpc to apply, if all the dark matter originates from particle decay. Since 
(AfWn > Age ~ 10 Mpc 3> 1 Mpc the bound (21) from non-dominance with b = 1 is 
naive. The actual lower bound on 5 lies closer to the hot dark matter (HDM) bound 
of Sec. 2.3.2 with 6 max < 1, cp. [59]. The heavier daughter from a two-body decay 
producing the desired dark radiation cannot form the observed dark matter. The HDM 
bound is the tightest bound obtainable taking into account the impact of the scenario 
on structure formation. 



Minimal velocity today Using (5) we can compute today's velocity v° of the decay 
products. Either a daughter is still relativistic, p > m v ~ 1, or has become non- 
relativistic, p < m v ~ p/m. With (4) and (3) we obtain for the heavier daughter 



o r o (^ + i) 2 -i (giy 

V ^ 6) = T d 2 5 + 1 [ajj ' ^ 

where <fi denotes the variance of the probability distribution derived in Appendix A. 
It is 4> ~ 0.215 if the decay occurs during radiation domination and ~ 0.376 if the 
decay occurs during matter domination. Bounds on today's dark matter velocity as- 
sume a Fermi-Dirac distribution, while the decay products are distributed as derived 
in Appendix A. While a transfer function is known [60], both distributions can be con- 
sidered as having an equivalent effect on structure formation as long as they have an 
identical second moment of the distribution (= variance), i.e. an identical root-mean- 
square velocity [61]. The factor cj) in (42) takes into account the exponential decay law 
in comparison to the sudden decay approximation. The minimal velocity of the heavier 
daughter today is obtained by inserting the lower bounds on 5 from Sec. 2.2 into (42). 
We find 

o km <f>(r) Ar /(U246\ ( <£,\% 

{V 2 )min ^ 3.3 — — — AA^cff - — -2 — • 43 

s 4>{t < t eq ) \Sldmh V \9* S J 

As the minimal free-streaming scale for r < t eq it does depend on the the amount of 
dark radiation only. Referring to [61] to suppress scales of the size of a galaxy cluster 
v° > 1 km/s. In this case bounds at least similar to the HDM constraints seem to apply. 
For decays later than £bbn we found a minimal velocity roughly three times larger and 
independent of the time of decay. If gt ~ 228.75 > g° s the velocity is significantly 
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lowered but still too large to be compatible with Lyman-a limits. This confirms that 
the non-dominance bound is naive and indeed HDM constraints apply in any case. 



2.3.2 Hot dark matter constraint and an opportunity 

If a particle is relativistic at matter-radiation equality and became non-relativistic in 
the meantime, it represents hot dark matter (HDM). In particular, observations of the 
structure in the Universe constrain the amount of hot dark matter, which counts to O m 
but does not form structures below its large free streaming scale. Observations are used 
to derive upper bounds on the sum of the masses of the SM neutrinos Y^ m u- Bounds 
vary between 0.4 eV for the minimal model including large scale structure data and 
2.6 eV if only CMB data is used and more parameters are included [62]. We re- write 
such bounds as constraints on the HDM fraction ^hdm/^dm < &max- It is 



S]d m V' ' / max p c m 



0.098 



(Tm u ) /eV, (44) 

V — / max 



where Odm — 0.21 has been inserted and n° = j^N^n® is today's number density of one 
neutrino species, if n® denotes today's number density of CMB photons. We note that we 
have not taken into account the differing phase-space distributions of neutrino hot dark 
matter, described by a Fermi-Dirac distribution, and decay products with a phase-space 
distribution given by exponential decay in an expanding universe, see Appendix A. There 
are differences in the damping tails between decay produced dark matter and warm dark 
matter [60]. Nevertheless, given current measurement uncertainties and considering that 
constraints on HDM arise from observations on larger scales we assume the HDM bounds 
to apply unaltered. 

As hot dark matter the corresponding decay product were non-relativistic today. 
Then its energy density p2 = m^Y^-so, when so denotes today's entropy density of the 
Universe and Y<i = nil s the particle yield. Its yield is the same as the yield of the 
decaying particle if it did not decay, Yi = Y. Thus its energy density as hot dark matter 
today is just suppressed by the mass ratio of the two particles. It is 

2 m 5+1 J 

where Qh 2 is given by (19). Any constraint on the amount of hot dark matter simply 
yields an implicit lower bound on S as in (25) with only ft^m — > ^hdm- ^ * s ''max < 1 
in (21) set by (44). For the bound depicted in Fig. 2 we choose a tight constraint, 
i.e. Yl m f < 0-44 eV at 95% CL [63] corresponding to 6 max = 0.043, because this 
yields a strong lower bound on 5. The temperature, when the heavier daughter becomes 
non-relativistic, corresponding to the obtained lower bound on 5 from hot dark matter 
constraints for r < (tJfOmin, i-e. T^id = 5^™), is obtained using (23) with 6 max given 
by (44). We see that T 2 nr ((5 = ffi fo) is lower than T eq , which is to be expected. The 
bound appears self-consistent without further assumptions. 
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If AiV e ff > is observed, while observations exclude the heavier daughter becoming 
non-relativistic before a certain time, the actual lower bound on 5 would stem from the 
requirement of the heavier daughter to become non-relativistic after this late time, for 
example, the time of photon decoupling. This possible constraint is depicted in Fig. 2 
as well. 

We point out an interesting fact that leads to an opportunity. The neutrino mass 
scale sets two in principle independent quantities. These are: i) the neutrino energy 
density after they became non-relativistic, 0, u ~ ^m v n v , and ii) the time neutrinos 
become non-relativistic. The CMB alone is not very much affected, if the neutrinos are 
still relativistic at the time of photon decoupling, cf. Sec. 6.1 in [64]. As the thermal origin 
and thus the neutrino temperature is understood, this restricts the CMB sensitivity to 
masses m v > 0.6 eV. The origin of HDM from a particle decay is very different. At 
what time the daughter becomes non-relativistic (8) depends on its mass only indirectly. 
The dependence is on the mass hierarchy to its mother as given by (45). Interestingly, a 
heavier daughter with an energy density above the tight HDM boundary from a decay 
with r <C t eq becomes non-relativistic at a time before photons finally decouple, < 
i^ pl . In this case we expect the CMB to be sensitive to the heavier daughter acting 
as HDM. Such an amount of relativistic energy becoming non-relativistic during CMB 
times should leave observable consequences in the CMB even though the observational 
situation is not clear anymore. For sure, observations on smaller scales are affected. 

By considering galaxy cluster data corresponding to a 10 Mpc scale iV eff = 3.91±0.42 
and J2 rn f = 0-34 ± 0.17 eV (both 95% CL) have been measured in [48] assuming free 
iV e ff and m u- This is even not at 2-a non-zero, but the maximum likelihood constraint 
is peaked away from zero. Improvements to this analysis are said to be already approved. 
The obtained mean value might well be the first hint of the neutrino mass scale, while 
the largest mass-squared splitting from neutrino data is y/\ Am^j ~ 0.050 eV [65]. In 
Fig. 2 the required 5 to have AN c s = 0.86 from the lighter daughter and fihdm = 0.004 
corresponding to Yl m u = 0.34 eV from the heavier daughter corresponds exactly to 
the upper edge of the thick solid black curve. Planck data combined with LSST or 
JDEM can constrain Y2, m v < 0.04 eV [66, 67]. Experiments like KATRIN [68] or those 
seeking the neutrinoless double beta decay [69, 70] will measure neutrino masses in the 
laboratory. We point out an interesting possibility given theses future sensitivities. If 
laboratory experiments measure m„ smaller than cosmological probes, this mismatch can 
be explained by a cosmological particle decay. The hot dark matter contribution in the 
Universe could have originated from particle decay. As we have shown the hot dark mat- 
ter and dark radiation can originate from the same decay. Since the connection between 
observed mass scale and the time, when the HDM becomes non-relativistic, is different 
than for an additional relativistic species, these two cases should be distinguishable by 
future cosmological observations. 
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Direct B™ x (t) from BBN constraints 
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Figure 4: Upper bound from BBN on the direct branching ratio of the decaying particle 
into electromagnetically interacting particles as function of its lifetime r. The thick solid 
curve represents the weakest bound obtainable. We consider larger branching ratios as 
excluded. The thin solid curve represents the strongest bound obtainable. It shows how 
strong the actual bound might become depending on the actual value of AN c s and the 
energy density of the decaying particle after t eq . Times are highlighted as in Fig. 2. 

3 Branching ratio constraints and opportunities 

In this section we derive constraints on several branching ratios of the decaying particle. 
Since the energy density of the decaying particle (19) is practically fixed for any signif- 
icant value AiV e ff > 0, these bounds are general and, in particular, independent of the 
particle physics model. 

Bounds from BBN We adopt the bounds from big bang nucleosynthesis (BBN) con- 
straints determined in [71]. There the direct decay of some hypothetical massive particle 
into pairs of Standard Model particles, X — > SM + SM, is assumed. The product 
f2 x Shad/em is bounded, when B denotes the branching ratio into either hadronically 
(had) or electromagnetically (em) interacting pairs of Standard Model particles. Elec- 
tromagnetic primaries are, for example, photons, 77, and electron-positron pairs, e + e~ . 
Hadronic primaries are quarks, qq, and gluons, gg. Hadronic decays generically lead to 
the injection of electromagnetically interacting primaries, for example, due to neutral 
pions decaying into two photons. The upper bound on corresponding branching ratios 
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is simply found as 

Qmax 

nmax _ had/em ,.„■. 
-°had/em ~ q J l 4D J 

where ^^^/ em is extracted from [71] and f2 is given by (19). We will argue that the 
bounds are indeed independent of the kind of coupling and the existence of only one 
dark decay mode. For example, they do not depend on 5, because 5 is in any case 
bounded from below by some cosmological requirement. So they might depend on this 
cosmological bound but not on the particle physics parameters. 

Bounds from BBN depend on the cosmic time. The desired energy density of the 
decaying particle depends on the time of decay as well. Therefore, the resulting bounds 
are characteristic functions of the lifetime of the decaying particle. They are unique 
for the production of cosmic dark radiation. Upper bounds from BBN on the direct 
branching ratio into electromagnetically interacting particles are depicted in Fig. 4. Since 
they arise from the destruction of formerly build nuclei by, for example, photofission, 
they become effective at rather late times ~ 10 4 s. However, they become severe with the 
strongest bound, B™ x ~ 1CP 5 , around 10 8 s. Afterwards the bound becomes weaker as 
the energy density of the decaying particle decreases following (18). Bounds are provided 
in [71] up to 10 12 s. Towards later times we perform a trivial linear extrapolation. Such 
an extrapolation is crude, but we will see that CMB constraints are anyway stronger in 
this regime. Since no AiV e ff > is confirmed, we have to consider a range of possible 
values. A minimal value AN^ n = 0.52 might be set by the expected 2-cr exclusion limit 
of Planck [19, 20]. As maximal value we take the 5-a exclusion of the combined analysis 
in [6]. It is AA^ ax = 5.265. These two different values lead to the spread of the two 
curves in Fig. 4 before the time of matter-radiation equality t eq . This is the same for 
the solid and dashed curves in Fig. 5. The different behaviour after matter-radiation 
equality in both figures stems from the following: Before t eq the conversion factor (17) is 
always nearly fj,, because 5 must be larger than <5 m j n given by (27). Here, we implemented 
the bound from the non-dominance requirement, because smaller 6 results at later times 
in a stronger strong bound. After i eq the maximal f2 is given by and 5 m ; n becomes 
smaller than one, which leads to the plateau of the strong bound after i eq . For the 
weakest bound obtainable the conversion factor (17) is /j, also after t eq , because S can 
also be arbitrarily large, if f2 ~ fi ra d instead. Thus the weak bound becomes weaker as 
Prad decreases relative to p mat due to the expansion of the Universe. 

Upper bounds from BBN on the direct branching ratio into hadronically interacting 
particles are depicted in Fig. 5. They arise from different processes at different times. 
Charged mesons and antinucleons affect relic abundances already at times as early as 
~ 10 _1 s. Therefore, the hadronic branching ratio is bounded, -BJJ^f ~ 10~ 4 -10~ 6 , 
already at such early times. The bound is strongest, ~ 10 -6 -10~ 9 , around 10 s 

and afterwards becomes weaker, because the energy density of the decaying particle 
decreases. In determining the hadronic constraints we encounter specific additional 
uncertainties: First, it has been shown in [71] that the hadronic BBN bounds depend 
not only on the time of decay, but also on the actual value of the hadronic branching 
ratio. We find that this dependence is too weak to be important for our purpose. 
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Direct B had (r) from BBN constraints 
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Figure 5: Upper bound from BBN on the direct branching ratio of the decaying particle 
into hadronically interacting particles as function of its lifetime r. Solid/dashed curves 
represent weakest/strongest bounds obtainable as in Fig. 4. Black curves apply to a 
mass of the decaying particle m = 1 TeV and grey (green) curves to m = 100 GeV. The 
thin dotted curves correspond in each case to a less conservative bound for the 6 Li/ 7 Li 
ratio. Within the enclosed area the cosmic lithium problems could be solved by the 
decay. Times are highlighted as in Fig. 2. 

Therefore, we assume a pure scaling of the bounds with the branching ratio between the 
extremal cases -Bhad = 1 and -Bhad = 0. Second, the bound on the hadronic branching 
ratio has also been shown to depend on the mass of the decaying particle. If the mass 
is varied while the energy density is kept fix, the number density varies accordingly. To 
first order different effects cancel out. Bounds for two different masses of the decaying 
particle are provided to indicate the remaining dependence. Black curves in Fig. 5 apply 
to a mass of the decaying particle m = 1 TeV and grey (green) curves to m = 100 GeV. 
Referring to Fig. 10 we note that it might be motivated by the production of dark 
radiation from particle decay to extend an analysis of BBN constraints towards smaller 
masses, m < 100 GeV, of the decaying particle. Third, the determinations of 6 Li and 
Li abundances are affected by uncertainties in the understanding of nuclei destruction 
processes in stars. The thin dotted curves in Fig. 5 correspond in each case to a less 
conservative bound for the 6 Li/ 7 Li ratio. Branching ratios above these curves but below 
the corresponding more conservative bounds should not be regarded as ruled out. 

In contrast, if they do not violate other bounds, they are a possible explanation for the 
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Figure 6: Upper bound from the CMB on the direct branching ratio of the decaying 
particle into photons and electron-positron pairs as function of its lifetime r. Thick/thin 
solid curves represent weakest/strongest bounds obtainable as in Fig. 4 from spectral 
distortions on the branching ratio into photons B^. In the same sense thick/thin dashed 
curves represent weakest/strongest bounds obtainable from the ionisation history on 
the branching ratio into electrons, positrons and/or photons. The dash-dotted curve 
indicates PIXIE's discovery reach, if AN c g = 1, Y p = 0.249 and g^ s = 2. Times are 
highlighted as in Fig. 2. 

relatively high 6 Li/H ratios observed in metal-poor halo stars providing the cosmic origin 
of 6 Li [45]. If the observationally inferred 7 Li/H ratio is solved by stellar depletion, both 
problems, known as the cosmic lithium problems could be solved by the same particle 
decay. At this point it is important to remind that the bounds in Fig. 4 apply equally 
to hadronic primaries, because they inject numerous electromagnetic primaries. These 
are also constrained in the following section. 

Bounds from CMB The emission of particles with Standard Model interactions leads 
to spectral distortions of the cosmic microwave background (CMB). As in the previous 
section, we exploit the fact that the energy density of the decaying particle is fixed by 
the amount of dark radiation. To determine an upper bound on the branching ratio into 
photons we update the analysis of [72] by taking into account the corrections pointed 
out in [73] and the current limits on deviations of the CMB from a thermal spectrum. 
These, obtained by COBE FIRAS, are |p| < 9 x 1(T 5 and y < 1.5 x 1CT 5 [74, 75]. We 
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find 

B, < 0.66 aw AN^g^g^U^* , (47) 

where p mSLX denotes the bound on p and i M ~ 6.91 x 10 6 s (1 - Yp/2) 4 / 5 the timescale of 
thermalization in the Universe at that epoch. It is set by the primordial helium abun- 
dance Y p and other cosmological parameters, where we inserted PDG mean values [51]. 
Following [73] we replace e^^ 5 4 — >■ 0.48(r/'^) 10 / 18 e L99( *' J// ' r ) 10/18 for times earlier than 
t/j. ■ Upper bounds on the direct branching ratio into photons from spectral distortions 
of the CMB are depicted as solid curves in Fig. 6. They are not effective at times earlier 
than 3 x 10 5 s, because injected photons thermalise safely not leaving any observable 
imprint. In the analytic approximation (47) this is explicated by the exponential factor. 
For t > tfj_ the weak bound becomes constant on a severe level, l?™ ax ~ 4 x 10~ 4 or 
2 x 10~ 5 . This is a qualitative difference to the BBN bounds. At a certain time around 
2 x 10 11 s they become stronger than the bounds in Fig. 4. Therefore, the extrapolation 
of the BBN bounds towards later times is not crucial for our purpose. The spread be- 
tween the strongest bounds obtainable (thin curves) and the weakest bounds obtainable 
(thick curves) is mainly due to the same reasons as for the bounds from BBN on the 
electromagnetic branching ratio depicted in Fig. 4. To consider the relatively large ob- 
servational uncertainty in the primordial 4 He abundance we take for the strong bound 
Y p = 0.267 and for the weak one Y p = 0.231, which corresponds to the 2-cr statistical and 
systematic PDG error range [51]. No change of slope arises at t eq from the change in the 
expansion law, because it turns out that in our parametrisation the bound function (47) 
has no proportionality to the time of decay. For times r < 4 VL^h 2 x 10 11 s elastic Comp- 
ton scattering establishes a Bose-Einstein spectrum with chemical potential p regardless 
of the details of the injection. For the considered case of massive particle decay, the 
number of injected photons is negligible relative to the number of photons in the back- 
ground. In any case the energy density of injected photons has to be small compared 
to the energy density of background photons. Then the induced chemical potential is 
proportional to the injected energy density, p cx B^p/p^. For times r > AU^h 2 x 10 11 s 
the spectrum can be described by the Compton y-parameter. The jump in the bound 
due to this change in the description of the spectrum with the corresponding constraints 
is easily identified in Fig. 6 around 10 10 s. It is B^p/p^ = Ay. 

Future CMB polarimeters such as PIXIE are proposed providing dramatically tighter 
constraints with projected detection levels of p ~ 5 x 10~ 8 and y ~ 10 -8 at 5-cr [76]. 
The dash-dotted curve indicates PIXIE's discovery reach, if AiV e fj = 1, Y p = 0.249 
and g^ s = 2. PIXIE could even identify the origin of a spectral distortion as particle 
decay [73, 77]. Comparing Fig. 6 and 9 with the corresponding hadronic bounds we 
see that due to the specific behaviour of the bounds large portions of parameter space 
probed by PIXIE are neither excluded by BBN nor by changing the ionisation history 
of the Universe. 

The emission of particles with Standard Model interactions may change the ionisa- 
tion history of the Universe, which can leave observable consequences in the CMB. We 
adopt the bounds on scenarios with late-decaying particles derived in [78] based upon 
WMAP7 limits. Upper bounds on the direct branching ratio into electron-positron pairs 
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Figure 7: As Fig. 4 but on the off-shell or non-tree-level branching ratio of the decaying 
particle as applying to the scenarios under consideration. Horizontal lines indicate i) the 
often quoted suppression of 10 -3 for off-shell processes and ii) the suppression from an 
electromagnetic loop and three-to- two-body kinematics. 



or photons from additional ionisation and heating observable in the CMB are depicted 
as dashed curves in Fig. 6. Xhey become effective at a much later time around 
and t eq , but also quickly stronger than the bounds from spectral distortions. They are 
strongest already around i^ pl and then become weaker towards later times as the energy 
density of the decaying particle decreases and a decay has less impact on observables. 
Nevertheless, they stay strong till today, in particular, stronger than bounds from spec- 
tral distortions at late times. In [78] for each lifetime it is scanned over the constraint 
for photons and electron-positron pairs for masses of the decaying particle ranging from 
2 keV to 12 TeV. The bound becomes a band with its width reflecting the variation 
in the bound between different decaying particle masses and decay products. To con- 
sider this uncertainty we take the strongest bound provided for our strong bound and 
the weakest bound provided for our weak bound, respectively. This approach is suffi- 
cient for our purpose, because the weakest bound obtainable from spectral distortions 
is already quite strong. Nevertheless, we would like to point out that this uncertainty 
-for example, 1.1 orders of magnitude at r = 10 14 s- could be reduced to a negligible 
level referring to [78]. Since we provide the injection spectrum and the variation with 
redshift, one should be able to do so using the provided grid of injection energies and 
redshifts. Referring to Fig. 10 we note that it might be motivated by the production of 
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Figure 8: As Fig. 5 but on the off-shell or non-tree-level branching ratio of the decaying 
particle as applying to the scenarios under consideration. Horizontal dashed lines indi- 
cate: i) the minimal hadronic branching ratio ~ 0.03 of a neutralino decay into gravitino 
calculated in [79] and ii) the minimal hadronic branching ratio ~ 10~ 4 of a sneutrino 
decay into gravitino for sneutrino masses > 200 GeV and 5 > 9 found in [80]. 



Bounds on considered scenarios In the previous paragraphs we determined cos- 
mological constraints on the direct decay, X — > SM + SM, into various decay prod- 
ucts with Standard Model interactions. As discussed these bounds are severe and any 
model must; satisfy them. For the scenarios considered in this work we assume that 
this is the case. In Sec. 2 we consider models without any tree-level decay mode that 
could give rise to a constrained branching ratio. This means that the constrained di- 
rect branching ratios are automatically zero. If, for example, photons are emitted on 
tree-level, the decay would have to occur before r < 10 4 s and so on. Thus all these 
constraints are fulfilled by construction. Dangerous summands in the Lagrangian are 
forbidden kinematics or by symmetries, which is relatively simple. This either holds 
true for the models under consideration in Sec. 4 or dangerous tree-level modes are 
sufficiently suppressed. We repeat that this often means branching ratios smaller than 
10 -4 . While tree-level branching ratios arise in a relatively simple way from the theory, 
higher order corrections contain further information about the theory. In any case de- 
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S vis (r) from CMB constraints in considered scenarios 
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Figure 9: As Fig. 6 but on the off-shell or non-tree-level branching ratio of the decaying 
particle as applying to the scenarios under consideration. 

cays emitting particles with SM interactions may proceed off-shell or via loop processes, 
X — > dark +?*—)• dark + SM + SM. Then the previously given bounds do not 
apply. Three-body final states might yield the leading contribution. Shortened, we have 
to consider how much energy is carried away invisibly in each decay. 

At early decay times, r <C t eq , the situation is very simple. Half of the energy will 
always be carried away by some dark decay product as we found already that both have 
to be emitted relativistically with equal momenta. The case that both decay products 
go off-shell or couple in a loop can safely be neglected, because the matrix element is 
additionally suppressed and leads to four-body final states. These are also suppressed 
kinematically. At later decay times, r > t eq , heavier decay products might be emitted 
non-relativistically. Actually, this is by construction the case for the strong bounds 
depicted in Figs. 4, 5 and 6. The virtuality of particles with p < m is small. It is very 
improbable that they decay off-shell or even couple in a loop. Therefore, it is a very good 
approximation to assume them to be safe. Only the decay products forming the dark 
radiation or -at least- being emitted relativistically are left to endanger observations 
in this case. Altogether, to determine upper bounds on the off-shell or non-tree-level 
branching ratio of the decaying particle we multiply (19) by ((5 + l) 2 — l)/(2(5 + l) 2 ). 
This is to consider the amount of energy which is always carried away invisibly. The 
second factor cancels the dependence on 5 in (19). Actually, this factor is universal in 
the sense of quantifying the relativistically emitted fraction of energy as we found that 



27 



m (t) with Y=F eq for DR production 
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Figure 10: Allowed masses m for a decaying particle with large equilibrium yield. Within 
the lower (mY = (mY) min , Y en (g* = gf M = 106.75)) and upper (mY = (mY) max , 
Yeq(g* = 228.75)) line a decaying particle, which left thermal equilibrium being rela- 
tivistic at some early time, might produce the desired amount of dark radiation depend- 
ing on and the type of particle. Horizontal lines indicate masses of the proton m p , 
the pion m^, the electron m e and photodissociation thresholds of 4 He and deuterium, 
respectively. Times are highlighted as in Fig. 2. 

5 is not free, but has to satisfy bounds. Applied to the scenarios of Sec. 4 this means 
Ba t = B\ = — ; ,n 9 , while the situation is odd, because it is restricted to an odd 
window of mass hierarchies, 1/4 < X2 < 1/2, see Sec. 4. The resulting bounds are 
depicted in Fig. 7, 8 and 9 analogously to Fig. 4, 5 and 6. Compared to those on direct 
branching ratios all bounds are reduced by 1/2. More importantly, there is no additional 
spread between strong and weak bound after i eq for off-shell and non-tree-level branching 
ratios. 

The obtained bounds are severe and have the power to exclude many particle physics 
possibilities. As demonstration we give two examples: i) From Fig. 4 a neutralino LOSP 
could decay during or even after BBN, r < 10 4 s, into a gravitino and a photon. The 
photon thermalises safely and the gravitino could act as dark radiation. However, as 
shown in Fig. 8 this is excluded for 4 x 10~ 2 s < r < i^ pl due to the minimal hadronic 
branching ratio of this decay. It arises from the off-shell decay of the massless photon 
into qq pairs enhanced by a logarithmic infrared divergence, see [79] for details. It is 
thus too large even for neutralino masses as small as some tens of GeV. ii) A sneutrino 
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LOSP can decay invisibly into a gravitino and a neutrino both possibly acting as dark 
radiation. The bounds in Fig. 4, 5 and 6 are evaded. However, the weak interactions of 
the sneutrino and the neutrino lead to a hadronic branching ratio with significant impact 
on BBN, if the sneutrino is sufficiently heavier than the electroweak gauge bosons. The 
minimal branching ratio of ~ 1CP 4 for sneutrino masses > 200 GeV and 5 > 9 found 
in [80] is depicted in Fig. 8. We see that the decay is excluded for 10 _1 s < r < 10 9 s. 

In Fig. 10 we show the band of allowed masses for a decaying particle, which entered 
thermal equilibrium and then decoupled being relativistic at some early time. Thus it 
has a relatively large yield. Larger yields can be reached in non-thermal production 
mechanisms. The figure serves to make an important point. We can see that for such 
yields the mass of the mother is restricted to be smaller than ~ 2 GeV for decays oc- 
curring during or after BBN. For masses of the mother smaller than the proton mass no 
nuclei can be emitted. Thus no BBN constraints from the injection of nuclei apply. For 
masses of the mother smaller than the pion mass no hadronically interacting particles 
can be emitted at all and so on. Emitted photons with energies below photodissociation 
thresholds do not destroy corresponding nuclei. A small enough mass of the decaying 
particle can thus circumvent all BBN constraints. It might be motivated by the produc- 
tion of dark radiation from particle decay to extend analyses of BBN constraints towards 
masses smaller than 100 GeV. There should be non-trivial effects, when the mass of the 
mother approaches the proton mass and so on. Regarding the ionisation history of the 
Universe the effect of photons emitted with energies smaller than one keV seem worth 
to be studied. 

4 Two dark decay modes 

In this section we study the origin of dark radiation from the two-body decay of a non- 
relativistic particle, which possesses two dark decay modes with corresponding branching 
ratios B\i summing up to roughly one, 

B(X -> 1 + 1) + B(X -»• 2 + 2) = B 1 + B 2 ~ 1 , (48) 

where 1 and 2 denote and label the two dark decay products. It is shown in Sec. 3 that 
the branching ratio into dark components is constrained to be very close to one at times 
later than £bbn- Considering the upper right corner of Fig. 1, there is no heavier/lighter 
daughter in each decay. In contrast, there are lighter daughters from one decay mode 
and heavier daughters from an additional decay mode. Compared to the cosmologies in 
Sec. 2 there is one additional parameter, the relative branching B1/B2. This allows for 
the dash-dotted curve in Fig. 1, where B2 is obviously much smaller than one. 

Well-motivated examples for such decays are i) saxion decays into two axinos with 
V 2 ~ m sax/( 4 /a)) where 77i sax denotes the saxion mass and f a the axion decay constant, 
ii) moduli decays into two gravitinos with y 2 ~ K 2 m|/(18Mpj), where denotes the 
modulus mass and k an effective coupling, iii) saxion decays into two axions with y 2 ~ 
x 2 m 2 SiX /(4:f 2 ), iv) flaton decays into two axions with y 2 ~ m 2 /(2f 2 ) [81], where mj 
denotes the flaton mass or v) moduli decays into bulk axions with y 2 ~ m|/(2M^j) [32, 
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33]. In i)+ii) a scalar decays into two fermions, while in iii)-v) a scalar decays into two 
scalars. Other combinations of spins are imagineable. Note that in theories beyond the 
Standard Model with a dark matter stabilising symmetry either interactions as thought 
about here or in Sec. 2 may be allowed for the same particle. Depending on the number 
of long-lived particles in the theory there may be particles of each kind in the spectrum. 
As one might infer from the given lists, this is the case, for example, in supergravity 
theories amended by the Peccei-Quinn mechanism. 

Even though all decay products may act as dark radiation, especially, the heavier 
ones do not need to do so at any time. In contrast, they may form the observed dark 
matter. Then two of three dark components originate from the same decay. This has 
been proposed first in [41], where the authors also take into account BBN constraints. 

4.1 Basics 

We introduce useful parameters, determine general properties and derive basic equations. 

Kinematics In each decay the two decay products have equal mass, so that the general 
expression for the initial momentum (1) of each emitted particle reduces to 

Pmi=Pi,2 = -(l-4x? i2 ) 5 (49) 

with X12 = mix/m denoting the mass ratio of the corresponding daughters with equal 
mass and the mother particle. This parameter is a useful measure of their mass hier- 
archy. We choose labels such that m\ < m<i- It is always mi, mi < m/2 otherwise 
the corresponding decay were kinematically forbidden. Non-relativistic emission can 
occur only in the odd mass window m/4 < m^ < m/2. This scenario might appear 
rather unattractive for model building from the naturalness point of view. The other 
way around, a theory predicting m<i in this range would interestingly always lead to this 
case. 

On \ and 1/x-t plane From the same general condition (7) we find the tempera- 
ture, when the corresponding daughters become non-relativistic, 

r 5 ^(^-4)-(|)' (50) 

with g™ = g* s (Ti T 2 ). Typical requirements on Tf\ might be that: i) the particle is 
still relativistic today, Tf r 2 < To, which is the best -or only- understood situation from 
the observational point of view, ii) the particle does not become non-relativistic during 
CMB times, rf r 2 < T 7 dpl or Jf T 2 > T(t™£), which would possibly leave some observable 
signature in the CMB, or iii) the particle does or does not act as radiation at matter- 
radiation equality, Tf \ < T eq or Tf\ > T eq . Such requirements turn into constraints on 
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Figure 11: Mass hierarchy l/x\ t 2 = m/mj^ given by (51) for different requirements on 
Tf\ ■ On the solid curve the decay products become non-relativistic at the time of photon 
decoupling. For smaller values of x particles become non-relativistic at later times. In 
the same sense the dashed curve corresponds to tf 2 = t™£, the thin dotted curve to 
tf T 2 = teq and for a decay product to be still relativistic today l/xi,2 would have to 
lie above the thick dotted curve. Within the horizontal dashed lines the corresponding 
decay product is emitted non-relativistically. Values of x > 1/2 do not exist. Times are 
highlighted as in Fig. 2. 



x\ 2 depending on the time of decay. From (50) it is straightforward to single out 

--#)?(£)'T- 

The constraints on x\p, from the requirements i)-iii) on Tf r 2 are shown in Fig. 11. 

Energy densities The desired amount of dark radiation determines the energy den- 
sity of the decaying particle. In principle, both decay channels can contribute to an 
increase in AN e g, which might be interesting on its own only if both add some observ- 
able contribution so for B\ ~ F>2- There is a plethora of possible cosmologies. For 
example, both can contribute dark radiation at BBN and/or photon decoupling. Heav- 
ier decay products could become non-relativistic in the meantime, while lighter ones are 
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still relativistic today, and so on. We shall focus on the case that the dark radiation at 
the time of observation is formed by the lighter daughters only, p dr = p\. Anyways the 
following holds analogously in and can be applied to other cases. 

Since B1 + B2 ~ 1, the number densities of the decay products at decay are given by 

ni i2 = 2nB 1>2 . (52) 

If the dark radiation is formed by the lighter daughter particles, 

Pdr = Pi = n\E\ , (53) 

where E\ ~ (pi), because they must be relativistic to act as radiation. At the decay p\ 
is given by (49) and n\ by (52). In this way the energy density of the decaying particle 
p = nm at the decay is set by the amount of dark radiation. The mean {pi)/pua = /i is 
given by the distribution in Appendix A. It is 

p\ dec = /i-^r 1 (1 - 4x?p Pdv (T d ) , (54) 

where pdr(T d ) has been found in (15). The previously defined conversion factor (17) 
would read in this case pB\(l — 4xf) 1 / 2 . 

After the heavier daughters became non-relativistic, their energy density, pi = 
U2E2 — n.2m,2, can be related to the energy density of the mother, if it would not 
have decayed, by making use of the definition of X2 and inserting (52). It is P2 = IB2X2P 
or today 

n 2 h 2 = 2B 2 x 2 nh 2 = 2B 2 x 2 ( ^) 4— . ( 55 ) 

VTd/ g* s pc 

where the energy density at decay is fixed by (54). Note that this is valid for times only, 
when the heavier daughters are non-relativistic. 

4.2 Non-dominance requirement and dark matter 

The non-dominance requirement, cp. (21), on the heavier daughters gives by (55) rise to 
an upper bound 



^_ £ , 626 x ,0-3, (1M) (|) 1 AN£ (1 _ ixl) i , (56) 



B. 



where we used further that B\^ are no independent parameters but linked by the def- 
inition of the scenario (48). In particular, B a /J3i = 1/{B^ 1 - 1). If then X2 or B2 are 
fixed by the particle physics model, the other of the two became fix making use of (56). 
For t > i eq the heavier daughters are either emitted relativistically as the lighter ones 
or their mass were restricted to the odd window m/4 < ni2 < m/2. For r < t eq the 
requirement simplifies, because: i) x\ is much smaller than one, cp. Fig. 11, and ii) B2 
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will be typically much smaller than one for the heavier daughter to make up a viable 
dark matter candidate, so that (B^ 1 — 1) _1 — B%. Then (56) reduces to 

A theory linking B2X2, such that this non-trivial constraint is fulfilled naturally, appears 
particularly attractive. 
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Figure 12: The required %2 = rri2/m depending on the time of decay to obtain a certain 
free-streaming scale A| s . The thick solid curve corresponds to the numerical solution 
of (58) for Af = 0.4 Mpc. The area with 0.2 Mpc < Af < 1 Mpc is found within the 
thin solid curves. Values above the upper thin curve are excluded if the heavier daughter 
is to form the observed dark matter. Values below the lower thin curve do not leave an 
observable imprint in the sky. The heavier daughters form cold dark matter in this area, 
if i?2 satisfies (56) at the boundary. The dashed curves provide the same information 
for the analytic estimate (60). Above the thick dash-dotted curve A| s > 10 Mpc. The 
thin dash-dotted curve corresponds to Af 3 = 100 Mpc and is given for better orientation. 
Overplotted as very thin grey curve is the simple analytic approximation (59). Barely 
visible for A| s = 4 Mpc it is slightly too small for 10 Mpc, while the fit for 100 Mpc is 
almost perfect again. We can read off that in a decay after 8 x 10 5 s the decay products 
have to be emitted non-relativistically to obtain Af s = 0.4 Mpc, because the required X2 
becomes smaller than 1/4. Other times are highlighted as in Fig. 2. 
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In principle after being emitted also the heavier decay products may act as dark 
radiation throughout the history of the Universe. However, they do not need to act as 
radiation at any time. At the boundary the heavier daughters form the observed dark 
matter, = ^dmi which energy density depends on the energy density of dark radiation 
as explicated in Sec. 2.3. We find that they can, indeed, form the observed dark matter. 
We will see how structure formation provides additional constraints on B2 and %i- To 
act as cold dark matter, for example, X2 must lie below the thin solid curve in Fig. 12. 
If this bound is violated, tighter constraints than (56) apply. For example, if they act 
as hot dark matter, we must replace fidm — > ^hdm- Since the free-streaming scale of 
the heavier daughters is set and can be adjusted by the involved couplings and masses, 
there are not only constraints but opportunities arising from the heavier daughters in 
structure formation. 



4.3 Solution to the missing satellites problem 



In simulated cold dark matter halos there is an overabundance of substructures with 
respect to the observed number of Milky Way satellites [46, 47] known as the missing 
satellites problem. The dispersion of structure on these small scales reduces the predicted 
number of galactic satellites. In particular, warm dark matter with a free-streaming scale 
Af s > 0.2 Mpc resolves the missing satellites problem [82, 83, 84]. We mentioned above 
the upper bound Af s < 1 Mpc from Lyman-a forest data. In [85] the authors considered 
the late decay of a massive particle with only one decay mode providing dark matter 
with the desired free-streaming scale. They found much smaller <5s than required by the 
consistent production of dark radiation, cf. Sec. 2. Their numerical results agree with 
those obtained from our formulae (38) and (41). 

We use the general result (37) for the free-streaming scale of a particle emitted with 
relativistic momentum. Inserting as below (37) with the only difference in T^ r , which is 
given, here, by (50), the free-streaming scale X 2 S of the heavier daughters as function of 
X2 and the time of decay is just (38) with the replacement ^-^j— - — > (x% 2 — 4)2. This 
can be read as implicit equation for x% arising from constraints on its free-streaming 
scale or to obtain a desired free-streaming scale to solve the missing satellites problem, 
for example. It is 



Af(x 2 ,r) 
0.4 Mpc 



nr\2\ 3 



+ 




(58) 



To obtain the thick solid curve in Fig. 12 we solved (58) numerically. Equally well, one 
can read off the required time of decay from Fig. 12, if x<i is set by the particle physics 
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model. Overplotted as very thin grey curve is the simple analytic approximation 

1.21 

z / u.4 ivipc \ 

X 2 ~ 
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10 7 s 
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(59) 
10 Mpc it 



For A| s = 0.4 Mpc it is barely visible except very close to 8 x 10 5 s. For Af 
is slightly too small, while for X 2 S = 100 Mpc the fit is again almost perfect. 

In order to find an analytic estimate we first simplify (58) by omitting rela- 
tivistic free-streaming and free-streaming after t cq . Exploiting logarithm rules we 
can see that one might neglect the factor (x 2 — 4)" 1 in the logarithm, if X2 3> 
(///2)(r/t cq ) 1 / 2 (^/(jr° r ) 1 / 4 ((7°g/^ s ) 1//3 , which holds for decays sufficiently earlier than 
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(60) 



where g™ = g* and g^ s = g^ s for X 2 S = 0.4 Mpc. As can be seen in Fig. 12 the result 
is systematically too large and the logarithmic dependence is misleading. On the other 
hand we find the dependence on r and a very weak dependence on g d . The shift with 
A| s is well reproduced. 

The corresponding time when the heavier daughters become non-relativistic is ~ 
3.4 x 10 6 s. So this event appears rather unobservable. It is not a surprise that we find 
a constant value, because in this case Af s = 0.4 Mpc is set by t™ as argued below (37). 
With the analytic result (60) we can verify the dependence of on g d . It is weaker 
than oc (<7*) 1/ ' 6 and mainly due to the change in the contribution from relativistic free- 
streaming. 

Inserting the desired X2 into the non-dominance requirement (56) at the boundary, 
so that the heavier daughters actually have to form the observed dark matter, yields the 
unique branching ratio B2 to produce the desired dark radiation and the dark matter of 
the Universe with the desired free-streaming scale at the same time from the same late 
decaying particle. We find a constant value, 

(61) 
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with an even weaker dependence on g* as i!J r . This is no longer a surprise, because with 
fix £| r a l so B2 becomes fixed in order for the heavier daughters to form the observed 
dark matter. Taking into account uncertainties 7.8 x 10~ 4 < B2 < 2.7 x 10~ 2 . In the 
course we had to discard the second solution to (58) as unphysical. To obtain an analytic 
estimate we insert (60) into (56) and single out 
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where we took 1 — Ax\ ~ 1 assuming sufficiently small x\ and, again, g^ T = g® and 
9*1 = 9*s- The logarithmic dependence is misleading, while the dependencies on AJ?, 
AiV e g and tldmh 2 in (61) find a reason. 

4.4 Hot dark matter opportunity 

In this section we point out the same hot dark matter opportunity as at the end of 
Sec. 2.3.2 but for the case of two dark decay modes. To be able to lead to a non- 
zero cosmological neutrino mass scale like ^ m v = 0.34 eV inferred from the cluster 
abundance as in [48] the free-streaming scale of the heavier daughters must be larger 
than the corresponding structure formation scale, A| s > A gc ~ 10 Mpc. As before we 
solved (58) numerically to obtain the required mass hierarchy X2 as depicted in Fig. 12 to 
have Af — 10 Mpc. The corresponding times for the heavier daughters to become non- 
relativistic are if(Af = 10 Mpc) ~ 2.4 x 10 9 s and f§*(Af = 100 Mpc) ~ 8.6 x 10 11 s. 
Thus only for very large free-streaming scales they become non-relativistic after t^^- 
Inserting the found X2 into (56) with the replacement f2dm — > ^hdm = b max (Y2 m^jfldm ~ 
0.033 f2d m we find at the boundary the corresponding branching ratio into the heavier 
daughters F>2 to obtain AiV e g = 0.86 and ^ m v = 0.34 eV from the same particle 
decay. In order to act as the dark radiation we assume the lighter daughters to be still 
relativistic today with the corresponding mass hierarchy x\ given by (51). We find a 
constant branching ratio 

B 2 ~ 10- 2 ( -^L_) ( AN, 1 ( 9*e¥.) . (63) 
VlO Mpcy V0.033y cfr V 0.12467 v ; 

Larger free-streaming scales are possible or mig ht be preferred. At Af = 100 Mpc the 
approximation (63) is significantly smaller than the true value ~ 0175. In any case, it 
is trivial to solve (56) with the appropriate replacements for any fixed value of X2- The 
allowed area of values X2(t) is found above the thick dash-dotted curve in Fig. 12. 

Altogether, we have shown that any desired amount of hot dark matter and dark 
radiation can originate from the same decay, while in the case of two dark decay modes 
the relative branching ratio allows to vary the time the hot dark matter becomes non- 
relativistic. This time determines the HDM free-streaming scale. It is relatively small, 
if they become non-relativistic before t™£ and larger, if they become non-relativistic 
after t™£. In the case of massive neutrinos and also in the case of HDM from particle 
decay with only one dark decay mode, this time is much later than t eq . This offers a 
possibility to differ between these cases in cosmological observations. 

5 Results and conclusions 

We studied particle decay as the origin of dark radiation. After elaborating general 
properties of such cosmologies we determined model-independent constraints on possi- 
ble underlying theories. Since the energy density of the decaying particle is fixed by the 
amount of dark radiation and thus by observations independent of an underlying parti- 
cle physics model, bounds on branching ratios and constraints on the mass hierarchies 
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between decaying particle and decay products depend in a unique way on the time of 
decay. 

If the decaying particle possesses only one dark decay mode, we find that the minimal 
free-streaming scale of its decay products is so large that hot dark matter constraints 
apply to their relic densities. Therefore, the heavier decay product in a particle decay 
producing the desired dark radiation cannot form the observed dark matter. So hot dark 
matter constraints determine the minimal mass hierarchy between decaying particle and 
the heavier decay product. This constraint is depicted in Fig. 2, where also cosmology- 
specific uncertainties are shown. The hot dark matter bound is the tightest bound 
obtainable taking into account the impact of the scenario on structure formation. On 
the other hand, a hot dark matter component in excess of the SM neutrinos could not 
only have originated from particle decay but also share its origin with the desired dark 
radiation. The heavier decay product can form a finite hot dark matter component, while 
the lighter one acts as dark radiation. In any case decay products become non-relativistic 
during or after CMB times. 

If the decaying particle possesses two dark decay modes, the free-streaming scales of 
the decay products are set and can be adjusted by the involved couplings and masses. 
Therefore, there are not only constraints, as shown in Fig. 11, but also additional oppor- 
tunities arising from the impact of the heavier decay products on structure formation. 
Depending on the time of decay we provide the unique mass hierarchy and relative 
branching to produce dark radiation and dark matter with any desired free-streaming 
scale from the same particle decay. The observed dark matter satisfying the cold dark 
matter paradigm may have originated from such a decay. Any finite hot dark matter 
contribution can be explained and possibly be distinguished in future cosmological ob- 
servations from other sources like SM neutrinos. In a different range of mass hierarchies 
and the corresponding branching ratios, the dark matter from particle decay solves the 
missing satellites problem. 

We determined general upper bounds on several branching ratios of the decaying 
particle into decay products with SM interactions. These are independent of the under- 
lying theory and, for example, independent of the number of dark decay modes. Direct 
decays are constrained as shown in Figs. 4, 5 and 6. Since these constraints are se- 
vere, we considered scenarios fulfilling them by construction. Dangerous terms in the 
Lagrangian are forbidden by kinematics or by symmetries. However, after taking into 
account the amount of energy that is always carried away invisibly we find that also 
off-shell and loop processes are severely constrained, cf. Figs. 7, 8 and 9. A certain finite 
branching ratio into hadronically interacting particles could solve the cosmic lithium 
problems. The emission of photons may enable future CMB polarimeters to detect and 
identify the desired cosmological particle decay. More robustly, the obtained bounds 
have the power to exclude many particle physics scenarios. As an example we show how 
decays of the lightest ordinary supersymmetric particle into an invisible particle like the 
gravitino are excluded as the origin of the desired dark radiation during and after BBN. 
We argue that particle decay as origin of dark radiation serves as a motivation to extend 
existing studies of constraints on cosmological particle decays towards smaller masses of 
the decaying particle. For example, BBN constraints do not seem to apply to a particle 
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that freezes out relativistically and produces the desired dark radiation in its necessarily 
late decay during or after BBN. 

Since we provide simple analytic formulae and figures pointing out uncertainties, our 
results can easily be adopted to constrain particle physics models and as a guideline 
for model building. Particle decay as the origin of dark radiation raises very specific 
requirements on any underlying theory. Most existing proposals assume implicitly that 
constraints are satisfied. Often decaying particle and decay products have only extremely 
weak interactions. Thus, they safely satisfy branching ratio constraints but cannot lead 
to additional observable consequences. 

We point out various opportunities of a cosmological particle decay serving as a 
motivation for further studies. There is a plethora of possible cosmologies to be explored. 
As close as some interplay between the different mysteries of our universe is, as attractive 
appears a theory that combines and intertwines them. 
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A Exponential decay in an expanding universe 

Any initial momentum from a decay is red-shifted by the expansion of the Universe. 
It is 



How the scale factor a oc t' i< - 1+u ''> grows with time depends on the equation of state of 
the Universe, p = up with u = 1/3 if it is radiation dominated and oj = in the case of 
matter domination. Earlier emitted particles experience a longer time of expansion and 
thus more red-shift than later emitted particles. Taking into account the exponential 
decay law this leads to a more involved momentum distribution than the monochromatic 
line obtained in the sudden decay approximation. 

Two-body decays have no intrinsic momentum distribution. In this case the momen- 
tum distribution function f(p,to) of an emitted particle is determined by the number of 
produced particles in a given time interval dt at time t. It is 



2 




(64) 



2 




/(p,*o) = -- r (66) 
r dp 
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x=p/p(r,*o) 



Figure 13: Normalised, time-invariant energy spectrum e(x = p/p(r, to); c) of relativistic 
particles from two-body decay (decaying particle at rest) in an expanding Universe 
with scale factor a oc t l ' c . The solid (dashed) line is obtained from (74) as described 
in the text with c = 2(3/2). The maximum of the energy spectra is highlighted at 
x = 1 ^ p = p(r,t ). 



and reversing (64) we can insert t(p,to), which gives dt/dp as well, also in N(t) arriving 
at 

f(p,t ) = cN p- c (T,t )e~(^) p c ~ l 

= cN p-^ ( ^-Y e-H^T , (67) 

where we abbreviated c = 3(1 + u)/2 > 0. With the assumption of a constant u we 
do not include the case of a particle that shortly dominates the energy density of the 
Universe at its decay during radiation domination. As argued in the text this situation 
is not expected to occur. 



Total energy Since dark radiation particles are relativistic, each particle's energy is 
given by its kinetic energy, so its momentum, E ~ p. The total energy in dark radiation 
is thus given by the integral 

Pini 

E dT (t ) = J pf(p,t )dp. (68) 
o 
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x=p/p(r,? ) 



Figure 14: Normalised, time-invariant probability distribution P(x;c) = cx c ~ 1 e~ xC 
of (76) for finding a relativistic particle from two-body decay (decaying particle at rest) 
within the infinitesimal momentum interval [p,p+Sp\. The solid (dashed) line is obtained 
for c = 2(3/2). Vertical lines highlight the maxima of the corresponding distributions, 
x ~ 0.707 (c = 2) and x ~ 0.481 (c = 3/2), as well as the mean of the distributions, 
/i = \/tt/2 (c = 2) and (i = T[2/3]2/3 (c = 3/2). The maximum of the corresponding 
energy spectra is at x = 1, cf. Fig. 13. 

Inserting (67) for c > this becomes 

£dr(io) = N p(T,t )- (t[-] - cT[l + -, ^] ) , (69) 

C \ C C T J 

where T[z] denotes the Euler gamma function and T[a, z] the incomplete gamma function. 
In the case of radiation domination (c = 2) (69) reduces to 



i . i 



2 t 



a 



E d r(t ) = N oP (T,t ) I ^-Erf[^J ] - e-^j , (70) 

where Erf[z] denotes the error function. 

Times of interest are later than the time of decay, so that we should investigate the 
limit to 3> T - I n this limit the general expression (69) yields 

lim E dr (t ) = N p(r,to)-T[-] . (71) 

to^>r c c 
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The difference by taking into account the exponential decay behaviour in contrast to the 
sudden decay approximation is thus found in the factor c _1 r[c -1 ]. For the two important 
cases this is 

[x = ~ 0.886 in radiation domination, c = 2, (72) 

and 

2 2 

\i = — r[— ] ~ 0.902 in matter domination, c = 3/2 . (73) 

3 3 

The difference in (71) between these cases is thus below two percent. However, compared 
to the sudden decay approximation it is twelve percent. 



Invariant energy spectrum The integral of the momentum distribution function 
f(p, to) over full parameter space yields the total number of particles gdr-^o- The integral 
of the energy spectrum e over full parameter space yields the total energy (68) and is 
thus just found as 

e(p,t )=pf(p,to). (74) 

The effect of red-shift is easily comprehensible. At times to ^ T it is an equal 
shift of all momenta towards smaller ones, which does not change the characteristic 
form developed roughly until to ~ 3r and to ~ 4r, respectively. Therefore, we give the 
normalised and time-invariant energy spectrum e(x) in Fig. 13 by defining x = p/p(r, to). 

We treat the non-relativistic decaying particle to be at rest. The correction from 
a non-zero kinetic energy is negligible, if the decaying particle has a sufficiently small 
momentum. This is the case for all considered scenarios in this work. 



Probability distribution Determining the probability distribution P(p, to) for find- 
ing a particle within an infinitesimal momentum interval [p,p + Sp] we note that by 
construction 

f(p,t ) 



i 



1= P(p,t )dp = -^r^dp. (75) 
Jo Jo iV o 

Thus we have also just found P(p, to)- Performing the same change of coordinates as for 

the energy spectrum, p — > x = p/p(r, to), we obtain 

roc 

P(x)dx= / cx c - l e~ xC dx . (76) 
Jo 

This time-invariant distribution is depicted in Fig. 14. Highlighted are the maximum of 
the corresponding energy spectrum at x = 1, which is independent of c, the maximum 
of the distribution itself at x ~ 0.707 (c = 2), x ~ 0.481 (c = 3/2) and the mean of the 
distribution 

l-OO 

H(P) = / xP(x)dx = c^Tlc- 1 } (77) 



to compare with (71). Obviously, P(x) is asymmetric. The variance 4>{P) also known 
as the second moment of P is found as 

, ■> f°° , ,9 , , 2, o f 0.215 for c = 2 

HP)= (x-^) 2 P(x)dx = r[l + -]-^^\ (78) 
Jo c 10.376 fore = 3/2 



41 



The difference in <j> is thus about 55% and compared to the sudden decay approximation 
like a factor of 5 in the case of radiation domination. 
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